Instability of flat disks with respect to the formation of twisted ribbons in smectic-A* 

monolayers 
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Smectic-A* monolayers self- assembled from aqueous solutions of chiral fd viruses and a polymer 
depletant can assume a variety of shapes such as flat disks and twisted ribbons. A first order 
phase transition from a flat disk to a ribbon occurs upon lowering the concentration of polymer 
depletant or the temperature. A theoretical model based on the de Gennes model for the smectic A 
phase, the Helfrich model of membrane elasticity and a simple edge energy has been previously used 
to calculate the disk-ribbon phase diagram. In this paper we apply this model to the nucleation 
process of ribbons. First, we study the "rippled disks" that have been observed as precursors of 
ribbons. Using a model shape proposed by Meyer which includes rippling in both the in-plane and 
out of plane directions, we study the energetics of the disks as functions of the edge energy modulus 
(a measure of the polymer concentration) and the mean curvature modulus k. We find that as 
the edge energy modulus is reduced the radial size of the ripples grows rapidly in agreement with 
experimental observations. For small enough k we find that the out of plane size of the ripples grows 
but its value saturates at a fraction of the twist penetration depth, too small to be experimentally 
observable. For large k the membrane remains flat though rippled in the radial direction. Such 
membranes do not have negative Gaussian curvature and thus will not likely spawn twisted ribbons. 
We also study the creation of twisted ribbons produced by stretching the edge of a flat membrane 
in a localized region. In experiments using a pair of optical traps it has been observed that once 
the membrane has been sufficiently stretched a ribbon forms on the stretched edge. We study this 
process theoretically using a free energy consisting of the Helfrich and edge energies alone. We add 
a small ribbon-like perturbation to the protrusion producd by stretching and determine whether 
it is energetically favorable as a function of the size of the protrusion. In qualitative agreement 
with experiment we find a nonzero value for the critical size of the protrusion needed to make a 
ribbon energetically favorable, though the value we find is an order of magnitude lower than the 
experimental value possibly due to our neglect of the director field. As in the case of the rippled 
disks we find that the mean curvature energy acts as a barrier between the disk and twisted ribbon 
structures. 

PACS numbers: 61.30.-v,61.30.Cz,64.70.M- 



INTRODUCTION 




Self- assembled Sm-^d* monolayers composed of chiral 
fd viruses in the presence of a polymer depletant have 
attracted much attention recently [l|, Q . An fd virus is 
approximately 880nm long and 6.6nm in diameter and 
its chirality can be tuned by temperature. These viruses 
behave like hard rods in the self-assemblies with a persis- 
tence length that is several times the virus length in the 
temperature range of the experiments [3]. Upon varying 
the concentration of the viruses or the polymer deple- 
tant, a variety of structures have been observed such as 
nematic tactoids, twisted ribbons and flat circular disks, 
with the latter two being Sm-^* monolayer membranes. 
A first order phase transition between large flat disks 
and twisted ribbons occurs upon lowering the concen- 
tration of the polymer depletant or lowering the tem- 
perature (thus increasing the strength of the chirality). 
During the transition from disks to ribbons, ripples start 
to form on the edge in the plane of the disk. The rip- 
ples grow in the in-plane radial direction as protrusions, 
making the membrane resemble the shape of a starfish 
as shown in Fig. [TJ To the accuracy of the experimental 
measurement (approximately 0.5 /am) the membrane re- 



FIG. 1: Creation of twisted ribbons as temperature is lowered 
(increasing chirality) [3]. At first ripples form at the edge 
of the membrane (leftmost figure). As the temperature is 
furthered lowered (going from left to right in the images) the 
protrusions transform into ribbons as shown. 



mains flat. As the transition proceeds, e.g., by continuing 
to lower the temperature, the protrusions transform into 
ribbons and consume nearly all of the material inside 
the disk. This transition is completely reversible upon 
increasing the temperature. It is also possible to trans- 
form a flat disk into a metastable twisted ribbon using 
optical traps attached at opposite edges to stretch the 
membrane [? ]. The membrane remains flat during the 
initial stage of stretching and then abruptly transforms 
into a twisted ribbon once the extension reaches a critical 
value as shown in Fig. ??. The resulting twisted ribbon 
relaxes back to a flat disk when the traps are removed. 

Theories describing the fiat disks and twisted ribbons 
have been constructed [l|, [H, using the de Gennes model 
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FIG. 2: Stretching of a membrane by optical traps [5] (only 
one trap, the bright spot, is shown). The membrane is in- 
creasingly stretched from one image to the next, from left 
to right. Once the membrane has been sufficiently stretched 
a twisted ribbon forms as can be seen in the two rightmost 
figures. 



for the Sm-A phase [8] generalized to include chirality, 
and in addition, in the case of the ribbons, the Helfrich 
model [il, [lo| for the surface bending energy. These the- 
ories use a simple form for the edge energy, the interac- 
tion of the rods at the edge with the polymer depletant, 
proportional to the edge length. A more realistic model 
incorporating surface tension and the "melting" of the 
smectic order at the edge has recently been developed 
prj . When applied to twisted ribbons the theory with 
the simple edge model yields good qualitative agreement 
with experimental measurements of the ribbon's pitch to 
width ratio providing that the Gaussian curvature mod- 
ulus appearing in the Helfrich energy is positive, in con- 
trast to the negative values typically measured in lipid 
monolayers or bilayers [l2|. By comparing the free en- 
ergy per unit area of a twisted ribbon with the corre- 
sponding energy of a large, flat membrane, a first-order 
phase transition between the two structures was found in 
agreement with experimental observation. However, the 
predicted value for the edge energy modulus was found to 
be an order of magnitude lower than that measured ex- 
perimentally, presumably due to the very simple nature 
of the edge energy model. 

The theory used to study the twisted ribbon is quite 
general and can be applied to a monolayer of any shape. 
The free energy F of the monolayer is given by: 



{fH^fn)dA^jd)dl 



where fn and fn are the Helfrich and de Gennes free en- 
ergy densities respectively and 7 is the edge energy mod- 
ulus ("line tension"). The Helfrich free energy density is 
given by 



fH = -K2Hf^kKG 



(2) 



where H and Kq are the mean and Gaussian curvature of 
the surface respectively, k is the mean curvature modulus 
(or "bending rigidity") and k is the Gaussian curvature 
modulus. We have assumed a zero spontaneous curvature 
because of the up-down symmetry of the fd viruses. 

The de Gennes free energy density in the one-elastic 
constant approximation with the assumption of perfect 
smectic order is given by [7]: 



-K 



(V • n)^ - 2qrv • (V x n) + (V x n)^ + 



(3) 



(1) 



where is the relative tilt angle of the director n with 
respect to the local surface normal, K is the single Frank 
elastic constant, q is magnitude of the spontaneous twist 
wave vector arising from molecular chirality and C is the 
tilt free energy modul us. The twist penetration depth is 
given by At = y^K/C. 

A monolayer of general shape can be modeled mathe- 
matically as a two-dimensional surface embedded in three 
dimensions. The surface is given by a position vector 
Y{ui^U2) parameterized by two coordinates, ui and U2. 
In terms of this position vector the total free energy 
Eq. ([1]) is found after some calculation to be |t7|]: 



^[{gikdjUk ^ giinkT^ji.)^ - (2nkLjk + cos6>)Yi]^ > v^g^^iig^^I2 



^ j s\v? 6^duidu2^k J ^ ^/g duidu2 ^ ^ J {g^^ LijY ^ duidu2 ^ ^ j) dl 



(4) 



where the indices j, /c = 1,2 and we sum over repeated 
indices. The tensors gij and Lij are the first and second 
fundamental forms of the surface, respectively; F^^- are 
the Christoffel symbols, g is the determinant of gij^ Cijk 



is the antisymmetric Levi-Civita tensor, and di = dm- 
The director field n is expressed in a local basis formed 
by Yi , Y2 and the layer normal N. More details on these 
quantities can be found in Refs. 0, fisl, fl^. 
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Henceforth we use dimensionless units measuring 
lengths in units of the twist penetration depth and 
energy in units of the Frank constant K. The twist pen- 
etration depth has been measured experimentahy with a 
value of 0.5/im[l]. The twist Frank constant in bulk fd 
systems has also been measured experimentally [HI with 
a concentration dependent value on the order of 10 ~^ 
dynes ~ lOO/c^T at room temperature. Thus, thermal 
fluctuations of the director are negligible, as noted al- 
ready in Refs. [l|,0l5 where the free energy was minimized 
and very good agreement was found between the pre- 
dicted director pattern and experimental measurements. 

In this paper we use this model to study the two mech- 
anisms for transforming a flat membrane into a twisted 
ribbon described above. In Sec. [Ill we study a model 
of the rippled disks observed when the transition to a 
twisted ribbon is driven by a decrease in the concentra- 
tion of polymer depletant which we assume is propor- 
tional to the line tension 7. We find that as 7 is de- 
creased the size of the ripples in both the radial and out 
of plane directions grow. In a narrow range of 7 there is 
a very rapid increase in the radial size of the protrusions 
while the height of the ripples remains nearly constant 
with a value too small to be observed experimentally. In 
Sec. inil we consider the transition to a twisted ribbon 
as induced by an external force, namely, a pair of op- 
tical traps which stretch the membrane (one trap holds 
the membrane in place while the second one pulls the 
edge of the membrane at the point of attachment). We 
consider a semi-infinite membrane and assume that the 
stretching produced by the second trap gives rise to a 
small protrusion at the edge of the membrane. We then 
add a ribbon-like perturbation to the tip of protrusion 
and determine whether the perturbation is energetically 
favorable. We find that for a large enough protrusion 
a ribbon-like perturbation is favorable. Due to the com- 
plexity of the shape we carry out this analysis for a model 
which neglects the director field and thus includes only 
the Helfrich and edge energies. We offer concluding re- 
marks in Sec. llVi 



II. INSTABILITY RELATED TO THE 
SPONTANEOUS PHASE TRANSITION 

We use a model [iG^ of the rippled disks where the rip- 
ples are described by sinusoidal waves both in the plane 
of the disk and out of the plane. While the experiments 
appear to observe fiat membranes this is only to an accu- 
racy of approximately 0.5 fim or one unit of dimension- 
less length. We allow for the possibility of out of plane 
fiuctuations in order to incorporate negative Gaussian 
curvature in the structure which we believe is necessary 
for the formation of twisted ribbons. We assume that the 
out-of-plane height decays exponentially into the interior 
of the disk so that the ripples are confined to the edge 
region. The disk lies in the x-y plane with its center at 
the origin. The height of the disk in the z direction is 



given by: 

/i(r, (j)) = Az sin(a^) exp[— 6(i?o — ^)]- (5) 

The radial coordinate i^(^) of the edge of the disk is given 
by: 

R{(j)) = Ro^Ar cos(a0 + V^). (6) 

Here r and (j) are polar coordinates in the plane, Az 
and Ar are the amplitudes of the ripples in the z- and 
radial directions respectively, Rq is the radius of the disk 
in the absence of ripples, a is the number of the ripples 
along the edge, b measures the decay of the height of the 
ripples into the interior of the disk and ip is the phase 
difference between the in-plane ripples and out-of-plane 
ripples. The radial coordinate r appearing in the height 
function /i(r, (j)) ranges from a cutoff ro 7^ up to R{(j)). 
For r < ro we assume that the disk is perfectly fiat. We 
introduce a nonzero cutoff ro because in the limit ro ^ 
the Gaussian curvature of the above shape diverges due 
to the factor sin(a0) in the height function. The value of 
ro is determined as follows. We calculate the Gaussian 
curvature ^^^(ro, Ar) for a ring of inner radius ro and 
outer radius ro + Ar with Ar small compared to Rq. 
We first choose ro of order Rq and then decrease the 
value of ro until we reach the point where Kciro^ Ar) 
begins to increase. This criterion determines the value 
of ro which is of order one for the shapes we consider 
here. Although the introduction of this limiting radius 
produces a discontinuity in the shape, the discontinuity is 
about 1% of Az and does not affect our physical results. 
Figure [3] shows an example of the rippled disk model. In 
Fig. m we show a contour plot of the Gaussian curvature 
of this shape. Note that the radial bulges are the regions 
of maximum negative value of the Gaussian curvature 
which we conjecture act as seed points for the growth of 
twisted ribbons. 

Although it is not clear from experiments whether 
these rippled disks are equilibrium structures, we assume 
they are and solve for the director field by minimizing 
the free energy, Eq. (|3|), computed for this shape. With 
an analytic form of the shape specified, the differential 
geometry quantities needed in Eq. dU can be computed 
explicitly. In principle we could then follow the approach 
of Ref. [7] and derive the corresponding Euler-Lagrange 
equations for the director field. However because of the 
complexity of the shape, the resulting equations are im- 
possible to solve explicitly even using numerical solvers. 
Instead we discretize the underlying x-y plane using a 
square lattice of grid size 0.05 and carry out an MC sim- 
ulation at low temperature (10~^ in dimensionless energy 
units with ks = I) varying the geometrical parameters 
Az and A^ of the shape. We note that our dimension- 
less energy unit, = 1, corresponds to approximately 
lOOksT at room temperature using the measured value 
of the twist elastic constant in fd solutions [l5| . 

The height of the membrane at each lattice site is given 
by Eq. ([5]) with Cartesian x-y coordinates converted to 
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FIG. 3: Shape of a rippled disk given by Eqs. ([5]) and (|6]) with 
= 0.5, Ar = 0.5, = 0.0, a = 5, 6 = 1.0 and = 5.0. 
(a) Viewed from above, (b) Viewed from a tilted angle. 



6F 




FIG. 4: Contour plot of the Gaussian curvature for the shape 
shown in Fig. [S] Note the maximum negative value at each 
outward bulge. 



polar coordinates. The shape of the membrane is fixed 
and not allowed to vary during the computation. Lattice 
sites with r > i^(^) (Eq. (|6|)) are excluded except those 
which have nearest neighbors lying within R{(j)) of the 
origin. These sites serve as "ghosts" which allow us to 
impose a boundary condition on the director field. As in 
Ref. we consider free boundary conditions [17], and 



the directors at the ghost sites are sampled in the MC 
simulation in the same fashion as the directors inside the 
disk. A finite difference algorithm was used to compute 
derivatives of the director field and the height function 
appearing in Eq. (j4]). In each MC cycle, the director 
at every lattice site in the disk and at the ghost sites 
was allowed to rotate with the magnitude of the rota- 
tion selected so that the total acceptance ratio was ap- 
proximately 50%. We started with a configuration where 
all the directors are pointing in the z direction and ran 
50000 MC cycles to reach equilibrium and then an ad- 
ditional 50000 MC cycles to collect data. Convergence 
was carefully checked. We have also verified that using 
a random initial condition for the director field produces 
the same equilibrium state. Unlike the usual method 
in a MC simulation where the free energy at each cycle 
is recorded and averaged, we recorded the orientations 
of directors in each cycle, computed an average orienta- 
tion for each director, and calculated the free energy for 
this averaged configuration. By doing so, we reduced the 
effects of thermal fluctuation and the conflguration ob- 
tained is close to the zero-temperature solution, i.e., the 
solution to the Euler-Lagrange equations. The validity of 
this method was checked by comparing the results of this 
method for a flat circular disk to the results obtained in 
Ref. |6| where the Euler-Lagrange equations were solved 
explicitly. 

We varied in Eq. (j5j) and found that for viruses 
with chirality q> {) the minimum of the free energy ap- 
pears at = as this value makes the handedness of the 
viruses and the edge of the membrane the same. Viruses 
with opposite chirality yield a minimum free energy with 
ijj = TT. Without any loss of generality we consider only 
the q > case. In order to study the instability to the 
formation of twisted ribbons, we vary the amplitude of 
the ripples, Az and A^, as the hue tension 7 (a measure 
of the concentration of polymer depletant) is reduced. 
We flx b at 1.0 in our dimensionless units. Physically, 
the length b is determined by the interplay between the 
deformation in the director fleld and the bending of the 
membrane so that one penetration depth is a reasonable 
choice. We flx the number of ripples a at either 3, 4 or 
5. We set the Gaussian curvature modulus ^ = 0.15 and 
the chirality q = 0.71 in accord with the earlier studies 
of flat disks and twisted ribbons [H, 0- Because the 
structure we are studying is not a minimal surface (i.e., 
one with zero mean curvature) we have allowed for both 
zero and nonzero values of the mean curvature modulus 
k as we discuss below. 

Fig. [5] shows the simulation result for the director fleld 
on a rippled disk with Az = 0.5, Ar = 0.5, a = 5, and q = 
0.71. Note that the directors are mostly perpendicular 
to the x-y plane at the center of the disk as expected. 

The results for the values of Az and A^ that minimize 
the energy as 7 is varied are shown in Table |I] for a disk 
with three ripples (a = 3) and k = 0. Note that the value 
of 7 found in Ref. [7] for the flrst-order phase boundary 
between disks and twisted ribbons is 0.276. We have 
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FIG. 5: Simulation result for the director field on a rippled 
disk with Az — 0.5, Ar — 0.5, a = 5, and q — 0.71. (a) 
Viewed from top. (b) Viewed from a tilted angle. 



TABLE I: Dependence of the height amplitude Az and the 
radial bulge amplitude Ar on the line tension 7 for a disk 
with three ripples. Note that value Ar — 2.2 in the bottom 
line is the largest value explored in our simulations yielding a 
lower bound on the value of Ar that minimizes the free energy. 
The values here are for the case k = 0,^ = 0.15. Within each 
range of 7 we have varied 7 with a step size of 0.001. 



7 


Az 




> 0.319 


0.0 


0.0 


0.314-0.318 


0.3 


0.3 


0.282 - 0.313 


0.3 


0.4 


0.274 - 0.281 


0.3 


0.6 


0.257 - 0.273 


0.3 


0.7 


0.246 - 0.256 


0.3 


1.0 


0.244 - 0.245 


0.3 


1.2 


0.243 


0.3 


1.3 


0.224 - 0.242 


0.3 


1.4 


0.216 - 0.223 


0.3 


1.7 


0.211 - 0.215 


0.3 


1.8 


0.200 - 0.210 


0.3 


2.0 


0.181-0.199 


0.3 


2.1 


< 0.180 


0.10 


> 2.2 



also studied disks with a = 4 and a = 5; our results are 
of the same order of magnitude as those shown in the 
Table for a = 3. The values for Az and Ar shown were 
obtained after exploring the range of values: < Az < 
0.5,0 < Ar < 2.2, a range selected on the basis of the 
order of magnitude of the experimental measurements. 
From Table H we see that when the line tension is large 
(well within the disk region of the phase diagram [7]), the 
free energy minimum corresponds to a disk where both 
the height and radial size of the ripples vanish. As 7 
is reduced the height and radial size of the ripples both 
grow. The height Az grows to a small value of order 0.3, 
which is well below the experimental resolution for out of 
plane fluctuations, and thus consistent with experimental 
observations of "flat" rippled disks. On the other hand, 
the radial size Ar grows more rapidly inside a narrow 
range of 7. Specifically, Ar quadruples from 0.3 to 1.2 in 
the range between 7 0.31 and 7 0.24, and continues 
to grow as 7 is reduced. These results are in qualitative 
accord with the experimental observations discussed in 
Sec. H namely, that ripples in the radial direction are 
observed and grow rapidly as the instability to twisted 
ribbons is approached. Recalling from Fig. |4] that ripples 
with nonzero Az have negative Gaussian curvature we 
conjecture that they are the seed points for the growth 
of twisted ribbons. 

The above results are for the case k = 0. Twisted rib- 
bons and flat disks are minimal surfaces with zero mean 
curvature; thus, the value of k is irrelevant to an anal- 
ysis of their energy. However, the present model of a 
rippled disk is not a minimal surface and thus we con- 
sider the effect of a nonzero value of k on our results. 
We find that for /c < 0.2 our results are qualitatively un- 
changed, however, the critical value of 7 where the rapid 
increase of Ar begins is lower. E.g., if A: = 0.2 we find 
that Ar reaches 2.2 (the largest value we have studied) 
at 7 ~ 0.175 compared to 7 ~ 0.180 when k = 0. For 
k > 0.2 we find that while the ripples grow in the radial 
direction as for the larger values of /c, the membrane re- 
mains flat, i.e., Az = and thus we presume that these 
ripples will not form twisted ribbons because of the ab- 
sence of seed points with nonzero Gaussian curvature. It 
appears then that the mean curvature energy acts as a 
barrier between the flat disk and ribbon states; a similar 
effect will be shown in the next section where we consider 
ribbon formation on a stretched flat membrane. 



III. INSTABILITY UPON STRETCHING 

In this section we model the instability of a flat mem- 
brane that is stretched using optical traps as described in 
Sec. [H Experimentally, this instability occurs within the 
region of the phase diagram where flat disks are the equi- 
librium shape. We consider large disks and use a semi- 
infinite model which has been shown [l| to accurately 
describe disks of the order of 10 jam in diameter. We as- 
sume that the optical tweezers produce a protrusion on 
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(a) (b) 
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FIG. 6: Model shapes of a membrane stretched by an optical 
trap before the ribbon- like perturbation is added. The colored 
region is the membrane. The Gaussian bump model is shown 
in (a) and the equilateral triangle model in (b). These figures 
correspond to zo = 2, A = 1 

the edge of the semi- infinite membrane and then deter- 
mine whether a ribbon-like perturbation attached to the 
tip of the protrusion is energetically favorable. Given the 
complexity of the analysis, we neglect the director field 
and consider only the Helfrich and edge energies. 

We assume the semi-infinite membrane lies in the x-z 
half-plane, z < 0. We consider two different shapes 
for the protrusion created at x = z = 0, as shown in 
Fig.ini The first. Fig. [6fa), is a Gaussian bump given by 
z = exp[— (x/A2:o)^] where A is a dimensionless con- 
stant setting the scale of the half width of the Gaussian 
bump. The second shape for the protrusion is an equi- 
lateral triangle with height zq as shown in Fig. [6fb). 

We denote the position vector of the stretched flat 
membrane (i.e., either of the shapes shown in Fig. [6]) 
by Yo(x,z). A small ribbon-like perturbation attached 
to the tip of the protrusion with a pitch axis along the z 
direction will produce a distortion of the surface in the 
direction of the surface normal N which lies parallel to 
the y axis. The position vector of the perturbed surface 
is then given by the following form: 

Y'=Yo + *(x,z)y (7) 

where ^{x^z) is the amplitude of the perturbation as- 
sumed to be small. Using the methods of Ref. Q we 
find the change 5Fh in the Helfrich free energy due to 
the perturbation: 

5Fh = 2kA ^kB^ 0(^^), (8) 

where 

^ = j {SHfdxdz = ^J (^,, + ^,,fdxdz (9) 

B = I {SKG)dxdz I (^,,^,, - ^l,)dxdz. (10) 

The subscripts on ^ denote partial derivatives taken with 
respect to the corresponding coordinate. The quantities 
dH and SKq denote the changes in the mean and Gaus- 
sian curvatures respectively due to the ribbon-like pertur- 
bation. In deriving Eq. [8]we have used the fact that the 



unperturbed shape specified by Yq{x,z) is a fiat mem- 
brane with H — Kg = 0. 

For a twisted ribbon with a pitch axis along the z di- 
rection and a pitch of magnitude 27r|6|, the y component 
of the ribbon's position vector is given by ^ = xtdin{z/b) 
[l3[. However, because the perturbation occurs only lo- 
cally around (0,zo) where the optical trap functions, we 
assume an exponential decay of the shape in the x direc- 
tion and in the z < direction (into the interior of the 
membrane). Thus, our final expression for the perturba- 
tion is given by: 



'^g(2; — 2;o)/Ai " 






b 


exp 


[-©1 



where Ai and A2 are parameters defining the rate of de- 
cay in the z and x directions respectively. The parameter 
b is kept large to ensure that ^ is small. Fig. [71 shows 
an example of the perturbation on a Gaussian bump pro- 
trusion. 

Substituting Eq. ([TTl) in Eqs. (j8j)-([l0]) yields the change 
in the Helfrich free energy, SFh- The change in the to- 
tal free energy (ignoring the director energy) is given by 
the sum of 6Fh and the change in the edge energy which 
is 7 times the change in the edge length. We evaluated 
the integrals in Eq. ([8]) and the change in edge length 
numerically. We assume that Ai = A2 = 1 using the 
following reasoning. If these lengths were small com- 
pared to the penetration depth, then the surface would 
be highly curved and there would be a large director en- 
ergy penalty. Conversely, if these lengths were large, then 
the perturbation would not be localized at the point on 
the membrane where the optical tweezer is active. We 
have verified that a small change in the values of Ai and 
A2 does not influence our qualitative results. Specifically, 
we have varied Ai from 0.1 to 1.5, and A2 from 0.6 to 4.0. 
We also assume that A = 1, as it is natural to expect 
that the half width of the Gaussian bump induced by 
the stretching is of the order of the extent of the stretch- 
ing, zq. We have verified that changing the value of A 
from 0.1 to 10 does not influence our qualitative results. 
Finally, we assume that b = 10^ and have verified that 
our results are to the accuracy of our calculation quanti- 
tatively identical as b is varied from 10^ to 10^. 

We set the Gaussian curvature modulus k = 0.15, the 
same value used in studying the phase diagram of fiat 
membranes and twisted ribbons [7] and considered sev- 
eral values of the mean curvature modulus k and line 
tension 7. Fig. [9] shows our results for the change SF in 
the total energy (Helfrich plus edge) as a function of the 
extent of stretching zq. In Fig. [9l^a), 7 = 0.3, which is 
in the flat disk phase of Ref. [7], and k = 0.0,0.05 and 
0.10 from bottom to top. In Fig. [D^b), k = 0.10 and 
7 = 0.3, 0.4, 0.5 and 0.6 from bottom to top (all in the 
flat disk phase). The dashed horizontal lines in both fig- 
ures mark the 5F = lines below which the perturbation 
is favored. For k ^ we see that the perturbation is not 
favored at the beginning of the stretching process, but 
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FIG. 7: Two views of an example of a ribbon- like perturbation 
Eq. on a flat semi- infinite membrane with a Gaussian 

bump protrusion on its edge. The parameters are Ai = A2 = 
A = 1.0, 6 = 10^ and zo = 1.778. Note the scale of the z axis 
in (a), consistent with the small value of ^. The view in (b) 
is along the y axis. Note that the projection of the structure 
on the x-z plane has the same shape as the Gaussian bump, 
Fig. EJa), consistent with Eq. (|7|. 



once a critical stretching zqc is achieved where SF = 0, 
the perturbation becomes energetically favorable. How- 
ever, the value zqc found in our model is one order of 
magnitude lower than the value of several penetration 
depths found in experiment. One possible source of the 
discrepancy is our neglect of the director energy. Re- 
call that the stretching experiments are performed in the 
region of the phase diagram where twisted ribbons are 
not energetically preferred in the absence of an external 
force. The theory Q used to study the transition from 
fiat membranes to twisted ribbons induced by lowering 
the concentration of polymer depletant showed that it 
is the chiral director energy that drives the transition. 
Thus, it seems reasonable to expect that including the 
director energy in the present case should increase the 
value of the critical stretching. This conclusion is further 
substantiated by noting that zqc increases with increasing 
line tension 7, i.e., increasing the concentration of poly- 
mer depletant, making the disks even more energetically 
favorable compared to twisted ribbons. 



FIG. 8: (Color online) (a) Integrated Gaussian curvature B 
(Eq. (dni), squares), integrated mean curvature squared A 
(Eq. (dni), triangles) and the change in the edge length 6L (cir- 
cles) as functions of the extent of stretching zo for a ribbon- 
like perturbation (Eq. (|lip ) on a Gaussian bump protrusion, 
(b) Same as (a) but restricted to the range < zq < 1. 

As can be seen from Fig. [9] (a), when k = the pertur- 
bation is favored at the very beginning of the stretching 
process, i.e., zqc = 0. Although disks and twisted rib- 
bons are both minimal surfaces with zero mean curvature 
where the value of k is irrelevant to the energetics, our 
present results suggest that the effect of k is to create a 
free energy barrier between disks and ribbons. This con- 
clusion is supported by the data shown in Fig. [8] where 
we plot the dependence of the integrated Gaussian curva- 
ture, B; the integrated mean curvature squared. A; and 
the change in the edge length, SL, as functions of the ex- 
tent of the stretching zq. Note that the integrated mean 
curvature squared term A remains nonzero (and positive) 
when the stretching is vanishingly small while the other 
terms (integrated Gaussian curvature and change in edge 
length) approach zero. Thus, the ribbon-like perturba- 
tion is not favored at the beginning of the stretching pro- 
cess because of the cost in mean curvature energy. While 
Fig. [9] (a) indicates that zqc grows as k is increased in 
value we find that there is a limit to this growth. For 
k ~ 0.6 we find zqc ~ 1.3. However, for larger values 
of k the ribbon-like perturbation is no longer favored no 
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FIG. 9: The change 6F in the total free energy (Helfrich plus 
edge energies) caused by a ribbon-like perturbation on the 
Gaussian bump as a function of the extent of stretching zo 
for different values of the mean curvature modulus k and the 
line tension 7. (a) 7 = 0.3, k = 0.0, 0.05 and 0.10 from bottom 
to top. (b) k = 0.10, 7 = 0.3,0.4,0.5 and 0.6 from bottom 
to top. The dashed horizontal line marks the line at 6F = 
below which the perturbation is energetically favored. 



FIG. 10: (Color online) (a) Integrated Gaussian curvature B 
(Eq. (dni), squares), integrated mean curvature squared A 
(Eq. (dni), triangles) and the change in the edge length 6L (cir- 
cles) as functions of the extent of stretching zo for a ribbon- 
like perturbation (Eq. (|11|) ) on a protrusion with shape of 
a equilateral triangle, (b) Same as (a) but restricted to the 
range < < 1- 



IV. CONCLUSION 



matter how large we make zqc- 

Turning to our second model of the protrusion, the 
equilateral triangle shown in Fig. [6fb), we display in 
Fig. [To] the dependence of A and 5L on the extent 
of stretching zq, similar to Fig. [8] for the Gaussian bump 
protrusion. The trends of the variations of these terms 
are similar to those in the Gaussian bump model. As 
in Fig. [Hl^b) we see from Fig. [TOT b) that the integrated 
mean curvature squared curvature is nonzero for = 0, 
supporting our contention that the mean curvature leads 
to an energy barrier between disks and ribbons. Similar 
to Fig. [9l the change in the total energy SF for the equi- 
lateral triangle protrusion is plotted in Fig. [11] for a num- 
ber of values of k and 7. For the triangle protrusion we 
find that 5F becomes negative when 2:0c = 0.724 ±0.001 
for k = 0.1 and 7 = 0.3 which is close to the value 
0.569 ± 0.001 found for the Gaussian bump. 



We have studied two kinds of instabilities of flat Sm- 
A* monolayers with respect to the formation of twisted 
ribbons. The first instability is related to the phase tran- 
sition from flat membranes to twisted ribbons which oc- 
curs when the concentration of depletant polymer (the 
edge energy modulus 7 in our theory) is lowered. We 
studied this transition using a model [16j of the rippled 
disks which are structures that have been experimentally 
observed as precursors of the transition to ribbons. Min- 
imizing the energy of the model shape we found that the 
size of the ripples remains very small when the line ten- 
sion is high (i.e., in the flat disk phase) while the size 
in the radial direction abruptly becomes large in a nar- 
row range of 7 when 7 is lowered. This result is con- 
sistent with experimental observations of the growth of 
the ripples and the rapid increase in the radial size Ar 
as compared to the imperceptible growth in the out of 
plane height A^. Ripples with nonzero A^ are regions of 
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FIG. 11: The change 6F in the total free energy (Helfrich and 
edge energies) caused by a ribbon-hke perturbation on the 
equilateral triangle protrusion as a function of the extent of 
stretching zo for different values of the mean curvature mod- 
ulus k and the line tension 7. (a) 7 = 0.3, k = 0.0, 0.05 and 
0.10 from bottom to top. (b) k = 0.10, 7 = 0.3,0.4,0.5 and 
0.6 from bottom to top. The dashed horizontal line marks the 
line at = below which the perturbation is energetically 
favored. 



relatively high negative Gaussian curvature and we spec- 
ulate that the rapid growth of with an accompanying 
nonzero value of is a signal of the instability to the 
formation of twisted ribbons. The second instability we 
studied occurs when a membrane is stretched using op- 
tical traps. We studied this phenomenon by considering 
a ribbon-like perturbation added to a protrusion created 
at the edge of a flat membrane by the optical trap. Our 
analysis of this instability was restricted to a free energy 
model which ignores the director field. Assuming that the 
mean curvature modulus k is nonzero we found that the 
ribbon-like perturbation is energetically favorable once 
the protrusion reaches a nonzero critical size zqc This 
result agrees qualitatively with experiment, however, our 
value of the critical stretching is an order of magnitude 
less than that observed experimentally, possibly due to 
our neglect of the director field. 

For both instabilities we have found that the mean cur- 
vature energy acts as a barrier to the creation of twisted 
ribbons. In the case of the rippled disk we found that if 
k > 0.2, the membrane remains flat and ripples are not 
energetically favorable. In the case of stretching the crit- 
ical value 2:0c grows with increasing k until k ^ 0.6 and 
zqc ~ 1.3. For larger values of k the ribbon-like pertur- 
bation is no longer favored no matter how large we make 
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